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OO ' I- INTRODUCTION 

As is known, dispersive shock is an oscillatory structure generated in wave systems after wave breaking of intensive 
. pulse under condition that the dispersive effects are much greater than the dissipative ones. In this sense, such shocks 
p I ' are dispersive counterparts of usual viscous shocks well known in dynamics of compressive viscous fluids. In the surface 
water waves physics, the dispersive shocks are known as tidal bores in rivers 0]- Besides this classical observation, 
the dispersive shocks have been also found in some other physical systems including plasma ^ and Bose-Einstein 
condensate 0,0]. 

In typical situations, a dispersive shock can be represented as a modulated nonlinear wave whose parameters change 
little in one wavelength and one period; hence the Whitham modulation theory 0,0 (see also 0]) can be apphed 
to its study. If one neglects dissipation, then dispersive shock is a non-stationary structure expanding with time 
so that at one its edge it can be represented as a soliton train and at the other edge as a linear wave propagating 
, with some group velocity into the unperturbed region. A corresponding Whitham theory of such shocks for the 
■ systems described by the Korteweg-de Vries (KdV) equation was developed by Gurevich and Pitaevskii in and 
I later it was extended to other equations such as the Kaup-Boussinesq system Benjamin-Ono (BO) equation 

and nonlinear Schrodinger equation 1141. This approach has found applications to water waves physics 
^§ I ^^'^ dynamics of Bose-Einstein condensate [l. Il6l| . 

. The Whitham method describes a long time evolution of the dispersive shock when many waves (crests) are 
^ ' generated. However, when the long-time evolution is considered, the small dissipation effects can become of crucial 
^ importance. In particular, they can stop self-similar expansion of the shock so that it tends to some stationary wave 
C structure which propagates as a whole with constant velocity. Correspondingly, the Whitham equations should be 
^ I modified to include the dissipation effects. For the first time it was done for the KdV-Burgers equation in [TtI IT^ 
by direct method which did not use applicability of the inverse scattering transform method to the KdV equation. 
More general approach based on the complete integrability of unperturbed wave equations was developed in (1^ and 
applied to the theory of bores described by the Kaup-Boussinesq-Burgers equation and the KdV equation with 
■ ■ ' Chezy friction and a bottom with a slope |2ll |. 

The method of Ref. 19 .1 can be applied in principle to any wave equation which is completely integrable in framework 
of the AKNS method [221 and any perturbation depending on the wave variables and their space derivatives. However, 
the important BO equation describing internal waves in stratified deep water includes the non-local dispersion term 
and therefore it cannot be considered by the method of Although the Whitham theory for the BO was discussed 
in [TH |23 | . its generalization to taking into account small dispersion effects has not been developed yet. The 
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aim of this paper is to develop the Whitham theory for the Benjamin-Ono-Burgers equation 

ut + uux + Hu^x ^ euxx, (1) 

where 

Hu[x) = -P.V. / ^^dy (2) 
TT J y~x 

is the Hilbert transform and the term in the right-hand side of ^ describes small friction with the viscosity parameter 
e. In the next Section we shall derive the Whitham equations which govern slow evolution of the nonlinear periodic 
wave due to its modulation and small friction and in Section III we shall apply this theory to the stationary bore 
(dispersive shock). 



II. WHITHAM THEORY FOR THE BENJAMIN-ONO-BURGERS EQUATION 

The unperturbed BO equation has a periodic solution 

u(x, t) — , V 0, 6 = kx ~ tot, (3) 

which depends on three constant parameters — wavenumber k, amplitude of oscillations A — {umax — Umin)/'2, and f3. 
This solution describes a nonlinear wave propagating with constant velocity 

V^'^ = IVA^+Ak'^ + p. (4) 
k 2 

In a modulated wave these three parameters become slow functions of the space and time coordinates and their 
evolution is governed by the Whitham equations which were obtained in for the general multi-phase solutions of 
the unperturbed BO equation by the method based on the complete integrability of the BO equation and in '13, 0| 
for the simplest one-phase solution by a direct Whitham method based on the use of the Hamilton principle 



with the Lagrangian 



dt J dxL{<j),(j),,^t)=0 (5) 



1 1 3 1 

L ^ -(l)t(j)x + -p:4ix + T:4'xH(f)xx, u^4>x. (6) 
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In this method the periodic solution is represented in the form 

(/.^V + '&W, (7) 

where 

ij) = I3x — ^t, 9 = kx — tot, (8) 

so that 

P^tpx, k = ex, uj^-Ot] (9) 

hence u = u(9,9x,0t,ipx): L = L{9,9x,9tT'ipxT^t)i and the averaging is taken over fast oscillations according to the 
rule 



2tt 



— I Ld9 (10) 



leading to the averaged Lagrangian which depends on the derivatives 9x,9t,'ipx,'4^t- The Euler-Lagrange equations for 
the corresponding averaged Hamilton principle 

d f dt f dxL{9x,9t,^x,^t)=0 (11) 
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yields with account of ® the Whitham equations in the form 

d_dL_d_dL_^ d_dL _ d_dL _ ^ 

dt dj dx 8(3 ' dt duj dx dk ' 

which should be complemented by the consistency conditions 

+ 7t = 0, fct + = 0. (13) 
After calculation of the integral H10() they reduce to the system of equations for the parameters l3,k,V ~ uo/k (see 

[Hill). 

Now our task is to generalize this procedure to the perturbed BO equation Instead of using the Lagrangian 
formulation with an additional field (see, for instance, Ref. 24]) we prefer to introduce a more simple approach 
which does not require introduction of new auxiliary fields. We propose to use directly the Hamilton principle in its 
infinitesimal form by noticing that Eq. can be written symbolically as 

dt j dx {5L + edp^rM} = 0, (14) 

where 5L is a variation of the Lagrangian ((HJ. Now we can transform ()14|l in the following way. First, we integrate 
the "dissipative" term by parts and use u = (f>x' 

J dt J dx{6L~euJu} ^0. (15) 

Second, we average the dissipative term as follows: 



dx 



dx{uxSu) = / dx[ux{ugSe + ug^59x + Ug^SOt + u^Jipx)] 

d 



UxU0 - -^{UxUgJ ~ — [UxUg^] 
OX Ot 



In the Whitham approximation with fast 6'-variable we have d/dx = 9xd/d9, d/dt = 9td/d9, where within the 
averaging interval the parameters 9x and 9t can be considered constant, the terms with ^-derivatives become equal to 
zero after averaging and, hence, we arrive at the expression 



dt j dx{uxSu) = J dt J dx{uxUg)59 = J dt J dxk{ul)S9. (16) 

Transformation of the term with Lagrangian in H15(l can be performed in a similar way and as a result we obtain the 
Whitham equations in the form 

d dL d dL „ d dL d dL 



dtdj dxd(3~^' dtduj axSfc ^^^^ 



which generalize Eqs. (|12|) to the BOB equation 
Simple calculation of averaged values gives [Tllll2| 



^-3^ -Hp-X+^ +6^ "2^^' (18) 



ek{ul)^2e(y~mV-PY~k']. (19) 

Their substitution into IjlTfl and use of (|13|l permit one to express 7 as 7 = /3^/2 and transform equations for the 
other parameters to the form 

A + PP. = 0, 

kt + {Vk)x = 0, (20) 
Vt + VVx + kkx = -e{V - mv - P? - fc']. 
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These are the Whitham equations for the physical parameters /3, k, V. 

Although equations (|20|l are simple enough for further investigations, they can be transformed to theoretically more 
attractive diagonal form by introduction of Riemann invariants a, 6, c according to definitions 

/3 = 2c, k^b-a, V ^ b + a, (c < a < 6), (21) 

so that we get the system 

at + 2aax = — 2e(a — c)(6 — c)(a + 6 — 2c), 

bt + 2bb.:^ = -2e(a - c)(6 - c)(a + 6 - 2c), (22) 
Ct + 2cCx = 0. 

In terms of Riemann invariants the periodic solution ^ takes the form 

u{x,t)^ + 2c, (23) 

a + 6-2c-2^(a-c)(5-c) cos 6* 

where 

9 = {b - a)x - {b^ - a^)t. (24) 

The parameters a, 5, c must satisfy the condition 2c < a + 6 to keep the solution non-singular. The amplitude of 
oscillations is expressed as 



A = A^Jia- c){b- c). (25) 

When for some concrete problem the solution of Eqs. H22|l is found and the functions a = a{x,t), b — b{x,t), c — 
c{x, t) are known, their substitution into (|23|l yields the modulated nonlinear wave for the problem under consideration. 
Mean value of u in this oscillatory region is equal to 

u = 2fc + /3 = 2(6 + c- a). (26) 

In the next Section we shall consider an example of such a problem. 



III. DISPERSIVE SHOCK (BORE) IN INTERNAL WAVES IN A DEEP FLUID 

As was mentioned in the Introduction, the dispersive shock is an oscillatory region joining two regions with different 
values of the wave amplitude u± which arise after wave breaking. In the simplest case of the Gurevich-Pitaevskii 
problem one can consider u± as two constants: u — s- u± as a; — > ±oo, respectively, and without loss of generality we 
can take u+ = and denote u_ = uq = const. Thus, we have to find the solution of the Whitham equations which 
corresponds to a modulated nonlinear wave satisfying the boundary conditions 

1^ 0, x ^ +00. ^ ' 

Naturally, the initial profile of the wave u = u{x,0) must also satisfy this condition and easy estimate shows that 
after wave breaking the waves are generated with wavelength L ^ I/uq. 
Now, we can distinguish two typical stages of evolution of the wave: 

1. The initial stage for t <C l/(eMo), when we have \ut\ ^ e|uj;^|, so that one can neglect a viscous term in ([^l; 

2. The asymptotic stage of large time t ^ l/(eMo), when the solution tends to the stationary solution determined 
by interplay of the dispersion and dissipation effects. 

Examples of the first stage have already been studied in 0, 0| . The simplest case of a step-like initial condition has 
been discussed in 11], and we shall reproduce some results here for convenience of future comparison with the second 
stage of asymptotically large time. 

Thus, we suppose that at the initial moment t = the region of transition from u = uq to u — is very narrow 
(i.e. much less than L ^ I/uq) so that the initial profile can be approximated by a step-like function. 



(28) 
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FIG. 1: Plots of the Riemann invariants &i t = 100 as functions of x in the case of decay of the step-like initial distribution 
without friction. 



As was noticed above, at the initial stage we can neglect the dissipation effects, so that arising dispersive shock is 
governed by the equations 

at + 2aax = 0, bt + 2bbx = 0, ct + 2ccx = 0. (29) 

After averaging over wavelength the initial- value problem for the Riemann invariants does not contain any parameters 
with dimension of length. Hence, the Riemann invariants can depend on the self-similar variable Q — x/t only and 
the Whitham equations reduce to 

a^(2a-C)=0, b^{2b~Q)^Q, c^{2c~Q)^Q. (30) 

Since at one edge of the oscillatory region we must have a — b (the soliton limit with fc = 0) and at the other edge 
a — c (the vanishing amplitude of oscillations limit A = 0), the only acceptable solution of these equations reads 
a = C/2 — x/2t, b — const, c = const. According to 128(1 . at the soliton limit x = x+ {a = b) the mean value u must 
match to the zero value at the right edge, u — 2c = 0, hence c = everywhere. At the opposite edge x = X- = the 
mean value must coincide with u = uq, u = 2b = uo, hence b = uo/2 everywhere. Thus, we arrive at the following 
solution of the Whitham equations (see Fig. 1) 

a — x/2t, < X < x+, 

b = uo/2, X < x-^-, (31) 
c = 0, X > X- = 

where 

X+ — Ugt. (32) 

Substitution of this solution into (|23|l gives the expression for u in the oscillatory region, 

( ,x {Upt - X)^ 

u[x,t) — — — -, (33) 

t{uot + X — 2y'uotx cos U) 

where 

= [2{uot - x)x - {ult^ - x^)]/{At). (34) 

The corresponding plot of the dispersive shock profile at fixed t is shown in Fig. 2. At the leading front we can see a 
soliton with the amplitude 

Us^W ^ Auq (35) 

which moves to the right with velocity V ^ uq. The trailing edge is located at x = and corresponds to a linear wave 
with vanishing amplitude and zero value of group velocity. Indeed, linearization of the BO equation with respect to 
small amplitude A in u Uq + Acos{kx — mi) leads to the dispersion relation 

to = uqU — . (36) 
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FIG. 2: Profile of the bore at t — 100 corresponding to tlie self-similar decay of the step-like initial distribution without friction. 



According to (|21|l wc have k = b = uq /2 at x = and hence 

= at a; = 0. 



dco 

dk 



k=uo/2 

Number of waves in the osciUatory region is equal to 

2 r^ot 



N 



2tt 



kdx = ^. 



(37) 



(38) 



The self-similar expansion of the oscillatory region holds as long as the viscosity effects can be neglected. However, 
these effects come into play at t ~ l/(eu§) and at i — > oo the shock profile tends to the stationary structure propagating 
with constant velocity. To find this structure, we look for the stationary solution of the Whitham equations H22|l so 
that the Riemann invariants are functions of (, ~ x ~ Vt only, where V — a + b ~ const. We assume that ^ = 
corresponds to the leading soliton front of the shock, where k = b — a = and u = 2(6 + c — a) = 0, which give at 
once that c = and a = 6=\^/2at^ = 0. Then the last equation H22|l gives c = identically and the rest equations 
reduce to a single equation 



(26- y)5{ = -2£V^6(V^- 6) 
which should be solved with the initial condition 



(39) 



b\i=o = (40) 

Elementary calculation with account of inequality 6 > a = y — 6, i.e. 6 > V/2, gives at once 6 = (y/2)[l + 
\/l — exp(2£T^^)]. At last, at ^ — oo we must have a = 0, b=V,u = uo = 2b= 2V, that is V = uq/2, and we 
arrive at the solution 



Uo 



1 - v/1 - exp(2eyO 



, Wo 



1 + v/l - exp(2eF0 



0, 



(41) 



where ^ — x ~ (uo/2)<. Plots of the Riemann invariants are shown in Fig. 3; they should be compared with Fig. 1. 
Substitution of Eq. 141() into Eq. (|23|l yields the profile of the shock. 



u{x, t) 



Mo[l - exp(£Mog)] 
1 — exp(euo^/2) cosO 



where 



Mo 



yjl - exp(euoOC, i 



X 1. 

2 



(42) 



(43) 



The corresponding plot is shown in Fig. 4. Again at the leading front we can see a soliton, but now it has the 



FIG. 3: Plots of the Riemann invariants as functions of ^ in the case of stationary bore solution with friction (e = 0.05). 



U(^) 




FIG. 4: Profile of the stationary bore described by the Bebjamin-Ono-Burgers equation with friction coefficient e = 0.05. 
Dashed line corresponds to the constant value of u at a: — > — oo. 

amplitude 

Us = 2uo (44) 

and propagates with velocity V = uq/2. Thus, the friction effects lead to decrease of the amplitude and velocity of 
the soliton compared with the non-stationary stage. However, it is important to notice that, on the contrary to an 
isolated soliton, the profile of the bore becomes asymptotically stationary. In this stationary solution the trailing edge 
is located at x = —oo. 



IV. CONCLUSION 

In this paper, we have discussed the structure of dispersive shock in the internal wave subject to small friction 
effects. It is shown that there exists a stationary profile so that non-stationary oscillating structure supported by a 
jump of the wave amplitudes at two spatial infinities tends asymptotically to this stationary profile. 

The Whitham method applied to this problem occurs very simple in the case of BOB equation and the corresponding 
equations can be solved in elementary and explicit form. 
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